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Abstract: We present next-to- leading order predictions for final states containing leptons 
produced through the decay of a Z boson in association with either a photon and a jet, or a 
pair of photons. The effect of photon radiation from the final state leptons is included and 
we also allow for contributions arising from fragmentation processes. Phenomenological 
studies are presented for the LHC in the case of final states containing charged leptons and 
in the case of neutrinos. We also use the procedure introduced by Stewart and Tackmann to 
provide a reliable estimate of the scale uncertainty inherent in our theoretical calculations 
of jet-binned cross sections. These computations have been implemented in the public 
code MCFM. 
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1. Introduction 

The production of vector bosons at hadron colliders provides a stringent testing ground 
for the Standard Model. The presence of a W or Z boson ensures that such processes 
are sufficiently hard to be reliably computed in perturbative QCD. The fact that the final 
states are produced via weak couplings means that it is a tough experimental challenge 
to observe these processes on top of the much larger QCD-dominated backgrounds. The 
search for processes involving multiple vector bosons also provides a crucial test of a key 
element of the Standard Model, namely interactions between the vector bosons themselves. 
Measuring these cross sections also allows one to place limits on additional self-couplings 
that may be induced by effective operators in many extensions of the Standard Model. 
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Theoretical predictions for the hadronic production of a Z boson and a photon have 
been significantly improved beyond the tree-level approximation. Next-to-leading order 
(NLO) QCD corrections [1-5] have been supplemented by higher-order gluon-initiated con- 
tributions [5-8] and NLO electroweak effects have also been included [9, 10]. These predic- 
tions, for the inclusive production of a Z7 final state, have been confronted with data from 
both the Tevatron [11, 12] and the LHC [13-15] experiments. However, particularly in the 
LHC environment, it can be advantageous to perform a less inclusive measurement in order 
to more cleanly isolate a potential signal. In particular it has become a common practice 
to separate an analysis into categories that are classified by the number of jets identified 
in the final state. Such an analysis is frequently referred to as a binned-analysis. This 
type of analysis can, for instance, be optimized to take advantage of the different expected 
background contributions in the various bins. Whilst this is a useful experimental tool it 
introduces new theoretical difficulties. This is due to the fact that the binning procedure 
introduces a new kinematic scale, namely the transverse momentum cut used to define a 
jet, upon which the theoretical prediction must depend. One normally expects this richer 
kinematic structure to lead to a less reliable perturbative calculation. 

At leading order (LO) the theoretical prediction for the Z7 cross section contains final 
states including only the Z and the photon. At NLO one includes virtual corrections to 
the LO topology and bremsstrahlung events corresponding to the emission of an additional 
parton. When the measurement is performed using exclusive jet-bins all of the NLO aspects 
of the calculation are confined to the 0-jet bin. Specifically the 0-jet bin includes the virtual 
corrections and the real emissions for which the parton is unresolved. The theory prediction 
for the remaining 1-jet bin is identical to performing a LO -Z^+jet calculation and ignoring 
the 0-jet bin altogether. In summary, the NLO calculation of the inclusive Z7 cross section 
provides a NLO prediction for the 0-jet bin, a LO prediction for the 1-jet bin and the 
prediction for bins with jet multiplicity of two or more is zero. 

The most natural way to improve the theory is to incorporate higher order perturbative 
corrections into the inclusive cross section prediction. This is a difficult task, although 
progress on this front may be expected in the relatively near future (see for example 
Ref. [16] for NNLO predictions for the diphoton process). In lieu of such a calculation one 
can instead focus on improving the theoretical predictions in the 1-jet bin. By performing 
a separate NLO calculation for the Zj+jet final state, the 1-jet bin can be predicted at the 
NLO level and a non-zero prediction is obtained in the 2-jet bin. This covers the kinematic 
range of the NNLO prediction for the Zj rate and has the same accuracy in the 1- and 
2-jet bins. Such a calculation does not, of course, include any higher order effects in the 
0-jet bin and therefore the accuracy of the Zj inclusive cross section is not improved. 

To this end, in this paper we present NLO corrections to the final state consisting of a 
Z boson, a photon and a jet. Although we will often refer to it in these terms ("^7+jet"), 
in fact we will actually compute the NLO corrections to the process, 

p + p-^l+7 + jet , (1.1) 

where the leptons are produced from either a Z boson or a virtual photon 7*. We will 
consider both charged and neutral leptonic decays i.e. t = e,v. When t is a charged lepton 
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Figure 1: Examples of leading order diagrams for the i^+jet process, for the cases of photon 
emission from the quark line ("g-type", left) and from the lepton line ("^-type", right). The particle 
labels correspond to momentum assignments that are all outgoing. 

we include the contributions in which the photon is radiated from the leptons. It is for 
this reason that the terminology "Z7+jet" is misleading, since it suggests the production 
of a Z boson, a jet and a photon, with the subsequent decay of the Z factorized from 
the process. As in the most recent study of the Z~f process [5], we include the effects of 
photon fragmentation in order to allow for photon isolation criteria that are currently used 
in experimental studies. 

As a by-product of this calculation, we shall also present results for the "Z77" process, 



where either, or both, of the photons may be radiated from a charged lepton. Many 
of the amplitudes relevant for this calculation can be obtained from the Z7+jet case by 
extracting the subleading-in-colour contributions. Although this process has previously 
been computed at NLO [17], we extend that treatment slightly by including fragmentation 
contributions. 

Our paper is structured as follows. In sections [2| and || we describe the analytic cal- 
culation of the Z^+jet and Z77 processes respectively. Section ||| discusses the issues of 
photon isolation and fragmentation. We present our results and some phenomenological 
examples at the LHC in Section || and summarize our findings in Section ||. Finally in 
Appendices [A| and || we present formulae for the helicity amplitudes that are used in our 
calculations. 

2. Calculation of Z^f+jet amplitudes 

In this section we present details of the analytic calculation of the NLO corrections to 
the Zj+jet process. Although NLO results exist in the literature for the W^7+jet and 
W77+jet final states [18, 19], the calculation of the ^7+jet process considered here is 
new. We include more detailed results in Appendix [A]. Note that we do not include higher- 
order finite contributions of the form, gg — > Zjg, that contribute at the level of a few 
percent at the LHC [20]. 



p + p — > l£ + j + 7 



(1.2) 
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Tree-level topologies associated with the hadronic production of a Z boson (with subse- 
quent leptonic decays), a photon and a jet are shown in Fig. |]. As illustrated in the figure, 
there are two separately gauge invariant classes of Feynman diagrams in the Zj+jet pro- 
cess. These can be classified according to whether the photon is: 

• emitted from the quark line, "(/-type" (Fig. [I], left), 

• emitted from the lepton line, "^-type" (Fig. [l], right). 

Clearly when the Z decays to neutrinos only the (/-type diagrams are present. Using the 
above nomenclature we can write the tree-level — > qqgjii amplitude in a form in which 
the explicit colour and electroweak charge structures are separated from the kinematics, 



A(°)(l 9 ,2 9 -,3 9 ,4 7 ,5,-,6 £ ) = 2V2e 3 [ g s T% 2 
x Q q {-Q q + v% R v{ )R V z (s 56 )) 4°)(l 9 ,2 g -,3 9 ,4 7 ,5,-,6 £ ) 

+Qi(-Q q + vl R vl R V z (U 56 ))Af\l q ,2g,3 g ,A y ,^Qi) • (2.1) 

Here the subscripts L and R refer to the handedness of the fermion that couples to the 
Z. The QED and QCD couplings are represented by e and g s respectively and Qi is the 
electric charge of particle i in units of e. The fermionic (quark/lepton) coupling to the Z 
boson, v£ is given by 

i -1 - 2Q f sm 2 9 w t 2Q^ sin 2 6>iy 



" u sin2(9n/ ' Vr sm2e w ' 

q _ ±1 - 2Q q sin 2 W q 2Q q sin 2 6 W 

Vl ~ S m26 w ' Vr ~ S m26 w ' 1 > 

where 6w is the Weinberg angle. The sign in v q L distinguishes between up (+) and down (— ) 
type quarks. The propagator factor, which is the ratio of the Z and photon propagators, 
is given by, 

VZ ^ = M*l T M ' 

s - M| + iT z M z 

where M z and T z are the mass and the width of the Z boson. We shall present expressions 
for the helicity amplitudes in standard spinor notation, with the spinor products defined 
as, 

(ij) = u^(pi)u + {pj), [ij] = u + {pi)u^(pj), (ij) [ji] = 2p, ■ pj . (2.4) 

A description of spinor helicity methods can be found in, for instance, Ref. [21]. The 
helicity amplitudes for the (/-type diagrams can be obtained from the primitive ampli- 
tudes for the e + e~ — > qqgg process presented in Ref. [22], once they are dressed with 
appropriately-changed color factors. At tree level this amounts to simply symmetrizing 
over the two gluon orderings in the partial amplitudes. For example, the tree-level ampli- 



tude A\ '(lq,2g ,3+,4+,5^,6/) can be obtained from Eq. (8.4) of Ref. [22] that reads, 

< ee (l+2r 3+4+) = - S^.J^l ■ (2.5) 
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Note that, following the original notation, we have suppressed the explicit dependence on 
the leptons in the amplitude definition on the left-hand side. Performing the symmetriza- 
tion one finds, 



■Aq 2 ? , 3,j~, Cjr , 5^ , 6^") — ^4g ree (l^, 2^ , 3^j~, 4+) + j4g ree (l^~, 2^ , 4+, 3jj~) 

(12) (25) 2 
\l3) (14) (23) (24) (56) ' 



(2.6) 



where the simplification results from use of the Schouten identity. Inspection of the above 
formula clearly confirms the correct QED structure. Compared to the original formula in 
Eq. ([T^), the QCD pole associated with the triple gluon vertex (34) has disappeared whilst 
new poles (13) and (24) have appeared. The new poles are associated with the fact that 
the photon is not colour-ordered. The remaining independent tree-level helicity amplitudes 
can be calculated using exactly the same prescription. The resulting expressions, together 
with the full details of the virtual and real radiation amplitudes, are given in Appendix [A|. 

The ^-type amplitudes cannot be extracted directly from the sub-leading colour pieces 
of the QCD amplitudes. This is obvious since the gluon can never be radiated from the 
leptons. However, as we shall explain in more detail below, we can take advantage of 
the fact that the QCD elements of the amplitude are all incorporated in the — > qqgll 
amplitude, and that the electroweak information (i.e. Z —> factorizes from the 

QCD part. Therefore we can use the e + e~ —> qqg amplitudes, given in the same notation 
in Ref. [22], and factor out the Z — > current, 

A q {l q ,2q,2, g ,h b Qi) = A^l q ,2 q ,3 g )J^(5 b 6 e ). (2.7) 

This extraction of the current from the amplitude is fairly straightforward. For example, 
consider the tree level helicity amplitude, 

2r,3+ 5^,6+) = -'(13/(23) (56) ' (2-8) 

In order to manipulate this into the desired form we first multiply both the numerator 
and denominator by [65] , in order to obtain a factor of s$q in the denominator that can be 
factored into the current J^. In the numerator, we use momentum conservation to write, 
(25) [56] = — <2|(1 + 3)|6] and obtain, 

w,^.*-jaz^ x (jm), ( , 9) 

where we have explicitly undone the Fierz identity. At this point we have identified the 
current, 

= (2.10) 

S 5 Q 



and the factorization in Eq. (2.7) is manifest with the identification 

I(1 + 3)|7m 
2(13) (23) 



^ 2p;) = i »l». (2 .n) 
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The £-type amplitudes can then be obtained by contracting the A^(l q ,2q,3 g ) piece with 
the Z —7- £ + £~j current, 

MU, 2?, 3 9 , 4 7 , 5,-, 6,) = ^(1,, 2 ? , 3 9 ) J^- 7 (4 7 , 5,-, 6,) , (2.12) 

where the currents for the Z boson decay including photon radiation from the leptons are, 

J ^ 7 ( 4 7< 5 7' 6 £ + ) = ^^J^^ef ' (2 ' 13) 
P* (A- k- «+\ [6|7 M I(4 + 5)|6] 

^ 7 (4 7 ,5,,6,)= , 456[45][46] • (2-14) 
Explicitly performing the calculation we obtain the following tree level ^-type amplitudes, 
•A-e ' '(lq "j 2= , 3j~, 4jr, 5j , 6j") = 2^ , 3j~) Jji „(4d", 5j , 6^~) 



(2|(1 + 3)| 7m |2) (5| 7 "|(4 + 6)|5) 
2 (13) (23) i 456 (45) (46) 
(25) 2 



9" 3+ 4" 5r 6 +N ) — A (\ + 2~ 3~h 7 M - (4~ 5- 6 + ) 

12 



(13) (23) (45) (46) ' 

-qi^qt ^g ) ° ft 

(2|(1 + 3)|6] ; 



t 456 (13) (23) [45] [46] 



where we have used momentum conservation in order to simplify the results. 

This procedure naturally extends both to other helicity amplitudes and to the NLO 
calculation. We have used the techniques described above to calculate the one-loop virtual 
amplitudes and real corrections to the Z7+jet process. In addition, the ^-type one-loop 
virtual amplitudes have been cross-checked with an independent calculation using analytic 
unitarity techniques [23-25], utilizing the Mathematica package S@M [26]. We present a 
more detailed breakdown of the calculation, as well as the full amplitudes, in Appendix [A|. 



3. Calculation of Z77 amplitudes 

We now consider the NLO corrections to the triboson process for which example 

leading order diagrams are shown in Fig. ||. NLO results exist in the literature both for 
this process [17] and VF77 [27], although this is the first time analytic results for this process 
have been written down. Comparing to the Z7+jet calculation presented in the previous 
section we observe that there is an increased number of distinct topologies related to the 
positioning of the two photons. The color structure of this process is trivial compared 
to the Z7+jet case, with only one such structure even at NLO. Using the same notation 
as the previous section we define the following sub amplitudes corresponding to the cases 
where; 

• both photons are emitted from the quark line (</</-type), 

• one photon is emitted from each of the quark and lepton lines (q^-type), 
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Figure 2: Examples of leading order diagrams in the Z77 process. The three diagrams correspond 
to emission of both photons from the quark line ( "gg-type" , left), one photon from each of the quark 
and lepton lines ("g£-type", centre) and both photons from the lepton line ("M-type" , right). The 
particle labels correspond to momentum assignments that are all outgoing. 



• both photons are emitted from the lepton line (££-type). 

Although there are nominally two qt topologies it is clear they are related to each other by 
the exchange of the two photon momenta. Explicitly, the decomposition of the tree level 
— > qqjj££ amplitude into our sub-amplitudes listed above is, 

A(°)(l 9 ,2 9 -,3 7 ,4 7 ,5 £ -,6 £ ) =AeH ili2 

x Q\ (-Q 9 + < i? <^z(s 56 ))4° ) ( 1 9> 2 <?>3 7 ,4 7 ,5,-,6,) 

+QiQ q (-Q q + vl R vl R V z (U 5e )) ^ ) (l g ,2 5 ,3 7 ,4 7 ,5 £ -,6 £ ) 



QiQq {-Q g + v q LjR vl iR r Z (t 356 )\ A$(l q , 2g, 4 7 , 3 7 , 5 



+Qj + < K 4,i?7M*3456)) ^2 ) (l g ,2 5 ,3 7 ,4 7 ,5^,i 



(3.1) 



The amplitudes for the qq-type diagrams are obtained from the subleading color contri- 
bution of the q-type diagrams in the -Z^+jet process, and similarly for the q£-type diagrams 
from the £-type Zj+jet. The amplitude for the ££-type diagrams may be obtained using 
the method that is described in Sec. ||, where we contract the qq — > Z/j* QCD current, 
Au(lq, 2q) with the Z/'j* — > £ + £~^^f electroweak current. However, in this paper we obtain 
the ££-type amplitudes simply from the qq-type amplitudes (for tree level) and from the 
g^-type (for tree level real emission) by crossing. The ££-type virtual amplitude is simply a 
vertex correction and is thus proportional to the ££-type tree level amplitude. The explicit 
relations for the amplitudes are presented in Appendix |B| 



4. Photon isolation and fragmentation 



Final states containing photons can provide useful tests of the Standard Model. Experi- 
mental analyses attempt to probe processes in which the photons directly participate in 
the hard scattering. However, such studies are complicated by the additional production 
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of photons through two mechanisms. Firstly, photons can be produced from the decays of 
unstable particles, for example 7r° — > 77. Since these photons are not produced directly in 
the hard scattering they are referred to as secondary photons. The second category occurs 
due to the fragmentation of QCD partons. Since the underlying production mechanism 
is a purely QCD process, these fragmentation photons are copiously produced at hadron 
colliders. 

Secondary and fragmentation photons are usually associated with significant amounts 
of hadronic activity. In order to reduce the effect of such photons, analyses typically require 
that the amount of hadronic energy in the vicinity of the photon is limited, 

#T ad < ^T ax with #0 = tJA^ 2 + At] 2 . (4.1) 

hade.Ro 



At the LHC, typical values for the parameters in Eq. (4.1) are a cone of size Rq ~ 0.4 and 



a limit on the maximum transverse energy, E^ ax ~ 5 GeV. This requirement is referred to 
as an isolation cut. 

Imposing such an isolation procedure raises an additional complication on the theo- 
retical side. Consider the production of a photon in association with a jet. At LO the 
matrix element is finite since the jet and photon cuts require that the two particles be well- 
separated in phase space. However, at NLO bremsstrahlung diagrams occur that involve a 
photon and two QCD partons. In addition to the usual QCD infrared singularities, which 
cancel in the combination with the virtual corrections, these diagrams contain a collinear 
singularity corresponding to photon emission from a quark. Since we require the photon to 
be resolved, this singularity has no virtual counterpart. Further, attempting to remove this 
singularity by requiring that no QCD radiation is present in the cone around the photon 
removes a region of phase space for soft gluon radiation, thus rendering the calculation 
infrared unsafe. 

One approach to avoid this problem is to isolate the photon in a different manner to 
that described above. An alternative isolation procedure, proposed by Frixione [28], allows 
arbitrarily soft radiation in the cone, whilst still eliminating the collinear fragmentation 
pole. The isolation criterion is, 

V E^ d 9(R - i? had)7 ) < e h E?J 1-008 * ) n for all R < R , (4.2) 
tad Vl-cosiW 

where n and are parameters of the algorithm. This prescription has the theoretical 
advantage that the fragmentation contributions do not have to be considered. On the 
other hand, this "smooth cone" isolation is difficult to apply experimentally. 

Alternatively, one can return to the definition of isolation that is used in the exper- 



iments, Eq. (4.1), and regularize the collinear fragmentation pole. This is achieved by 
absorbing the collinear splitting into the kernel of a fragmentation function that behaves 
in a manner analogous to initial state PDFs [29]. These fragmentation functions, which 
satisfy the DGLAP equation, must be extracted from data due to their non-perturbative 
nature. The NLO prediction then consists of the usual NLO diagrams plus a QCD LO 
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matrix element coupled to the relevant fragmentation function, 



a 



.7 

NLO 



= a 



.7 

NLO, direct 




(4.3) 



a 



In this equation a°j represents the production of a final state parton of species a in place 
of the final state photon, and D a ^(z, Mp) is the fragmentation function with p 7 = zp a . 
Note that the introduction of the fragmentation functions also requires the addition of 
a new scale, the fragmentation scale Mp. As a result the separation between the direct 
and fragmentation pieces is not unique and only the sum is theoretically well defined. In 
MCFM we have implemented the fragmentation sets of BFG [29] and GdRG [30] . 

5. Phenomenology 

We have included the processes described in sections || and |3] into the NLO parton level code 
MCFM which is available publicly [31]. This builds on the existing Wj and Zj processes 
already in the code. MCFM uses the dipole subtraction scheme formulated by Catani and 
Seymour [32] in order to isolate singularities in the virtual and real contributions. We use 
the following MCFM default electroweak (EW) parameters in our calculation, 



The remaining EW parameters are defined using the above as input parameters. In our 
calculations we use the CTEQ6L1 PDF set at LO, the CT10 PDF set at NLO [33] and Set 
II of BFG [29] for the fragmentation functions. 

5.1 Zj + jet at the LHC 

In this section we investigate the phenomenology of the production of Z7 and an associated 
jet at LHC operating energies. For typical experimental photon selection cuts (with pip > 
15 GeV), the e + e~j inclusive cross section is relatively large, about 1 pb at 7 TeV [5]. As a 
result, the process in which an additional hard jet is radiated can also be readily observed. 

We therefore begin by presenting ££j+jet cross sections for the LHC at a range of 
operating energies. We base our selection cuts on those used in the most recent ATLAS 
analysis [15]. In order to pass the selection cuts an event must satisfy the following criteria, 



Photon : |%| < 2.37, R h > 0.7 , 

pip > 15 GeV (low— pt), 

pip > 60 GeV (intermediate— pt) ; 

pip > 100 GeV (high-p T ). 
Leptons : m t + t - > 40 GeV, p e T > 25 GeV, |%| < 2.47, R tj > 0.6, 

Ep iss > 25 GeV . 

Jets : pip > 30 GeV, \yj\ < 4.4, k T algorithm with R = 0.4 . (5.1) 



M z = 91.1876 GeV . 
T z = 2.4952 GeV 
G F = 0.116639 x 10~ 4 GeV 



-2 



M w = 80.398 GeV , 
T w = 2.1054 GeV , 
, m t = 172.5 GeV . 
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Photon Cut 


Cross Section 


= 7 TeV 


= 8 TeV 


loW-pT 

04 > 15 GeV) 


LO 

e + e~7 + jet + X 

NLO 


1 ae+ll-6% 
ioo -10.0% 

9 cro+7.4% 
ZOZ -5.2% 


OOQ+9.9% 

zzo -9.4% 

q m +5.9% 
OUi -5.4% 


intermediate-pT 
(pi > 60 GeV) 


LO 

e + e~7 + jet + X 

NLO 


1 8 R+16.5% 
i0 ' D -13.3% 

24 6 +8 - 1% 

Z ^' U -6.8% 


nq n+15.4% 

zo.u_ 12 6% 

qn 7 +6.9% 
^•'-6.8% 


LO 

3(W)7 + jet + X 

NLO 


i S9 +16.8% 
iOZ -13.5% 


99Q +15.7% 
zzy -12.8% 

qi /?+9.2% 
01D_ 79% 


{pi > 100 GeV) 


LO 

e + e~7 + jet + X 

NLO 


fi 7n +19.0% 
°-' U -14.9% 

O q/,+6.6% 


co+17.9% 
S - 0Z _14.4% 

1 n 7 +6.0% 
lu -<_6 .1% 


LO 

3(z/p)7 + jet + V 

NLO 


65.4^-;^ 


cm fi +18.0% 

i i 9 +8.9% 
11Z -7.8% 



Table 1: LO and NLO cross sections (in femtobarns) for e + e _ 7+jet and 1^7+jet production 
at the LHC with centre-of-mass energy 7 and 8 TeV. Note that, for the neutrino case, we have 
summed over three flavours of neutrino. The quoted value represents the scale choice of Mz with 
the percentages indicating the shift in the central value upon varying this scale by a factor of two. 

In addition photons are isolated by requiring that E^ A G R < £f? ax with R = 0.3 
and £?p ax = 6 GeV. LO and NLO cross sections under the fiducial cuts listed above, for 
operating energies of 7 and 8 TeV, are shown in Table |]. The factorization, renormalization 
and fragmentation scales are set equal to one other, = I^R = Mp = Mz- Note that we 
have not provided predictions for the neutrino case with a low photon cut since it may not 
be possible to trigger on such events. First, we see that the effect of tightening the photon 
cuts from 15 to 60 GeV is to lose about an order of magnitude in the yield. Raising the 
cut even further, to 100 GeV reduces the cross sections by a further factor of three. In 
all cases the effect of the NLO corrections is to increase the cross section by a factor of 
about 1.3, independent of the operating energy. For the e + e~7+jet final state using the 
low photon pt cut, we extend the dependence on operating energy up to the LHC design 
target of 14 TeV in Fig. ||| As is clear from the figure, the JC-factor is constant across the 
entire foreseeable range of LHC operating energies. At 14 TeV the cross section is twice as 
large as at 7 TeV but the fractional uncertainty is approximately the same. 

In this section we have focused on NLO results for inclusive quantities. In the following 
section we will more closely follow the experimental setup by splitting the analysis into jet 
bins. 

5.2 Z7 and Z7+ jet : Exclusive predictions 

In this section we will investigate more exclusive quantities involving a final state consisting 
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Figure 3: Dependence of the LO (blue) and NLO (red) inclusive cross section of pp —> £ + £ 7+jet 
process on the renormalization/factorization scale, for a range of LHC operating energies 



of Z^f and a fixed number of jets. As we discussed in the introduction, previous versions of 
MCFM [5] are able to predict quantities accurate to NLO in the 0-jet bin, and to LO in the 
1-jet bin. Using the calculations presented in this paper we are able to extend the 1-jet bin 
to NLO accuracy. With the ATLAS results in hand [15] we will also take the opportunity 
to re-evaluate predictions for the 0-jet bin and reassess their theoretical uncertainties. 

The ATLAS paper [15] studied the production of Z7 final states, separating events 
based upon the number of reconstructed jets. Using this approach they were able to provide 
measurements of both the inclusive i.e. a{V^ + X) and exclusive i.e. a{V^f + jet) cross 
sections. Explicitly, the results for the £ + £~ , j final state (averaged over electron and muon 
channels) quoted in Ref. [15] are: 

Low— pt region : 

Inclusive : a ATLAS = 1.29 ± 0.05 ± 0.15 pb 
Exclusive 



a ATLAS = LQ5 ± 04 ± 0.12 pb 



a MCFM = 1.22 ± 0.05 pb 

MCFM 



a 



1.03 ± 0.04 pb 



Intermediate— pt region: 

a MCFM = 58±5fb 



Inclusive : a ATLAS = 68 ± 8 ± 5 fb , 
Exclusive : a ATLAS = 47 ± 7 ± 4 fb , 



a 



■™ = 40i3fb 



(5.2) 



The experimental measurement is given with the statistical (first) and systematic (second) 
errors separately. The MCFM prediction is taken from the parton level values quoted in 
the paper (the collaboration also presents a particle-level corrected result). We see that, 
although the low-p^ results are in agreement with the MCFM prediction, the intermediate- 
Pt measurements are somewhat higher, albeit with larger uncertainties. A further feature 
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Figure 4: A comparison of estimating scale uncertainties using the standard method of varying 
around a central scale by a factor of two (red) with the ST method (blue). The LHC operating 
energy is 7 TeV and the photon cut takes the low value, i.e. p7p > 15 GeV. The left hand plot shows 
the error estimation at NLO in the exclusive 0-jet bin. The plot on the right side represents the 
same quantity for the NLO 1-jet exclusive cross section. 



is also apparent: the relative theoretical uncertainty on the exclusive cross sections is com- 
parable, or smaller, than that on the inclusive prediction. Since the exclusive calculation 
introduces a new scale, corresponding to the transverse momentum of any additional jets 
that are effectively vetoed, one expects the calculation to have a richer structure and hence 
a larger theoretical uncertainty. The principle aim of this section is thus to investigate 
an improved method for estimating the scale uncertainty [34] in order to provide updated 
theoretical results for the 0-jet bin. In addition we will provide new theoretical predictions 
for the 1-jet bin that can be compared against future measurements. 

The naive method of estimating theoretical uncertainties by scale variation, which is 
of course only a crude estimate of missing higher order contributions, has been shown to 
be extremely dangerous in the presence of jet vetoes [34]. The authors or Ref. [34] propose 
that instead of using the scale variation of the exclusive cross section as a measure of the 
uncertainty one should use the following, 

A% = Al Zj + A 2 > Z7j . (5.3) 

In this equation Az-y represents the total uncertainty in the exclusive 0-jet bin and A>z-y 
and A>z-yj represent the uncertainties obtained from the inclusive calculation of Z'y and 
Z7+jet respectively. In this way the two perturbation series in as which make up the 
exclusive prediction are treated as uncorrelated. This ensures that no accidental cancel- 
lation between the scale-dependent coefficients in the perturbation series occurs. As is 
stressed in Ref. [34], if A>^ 7 is calculated at NLO then A>z-yj should be calculated at LO. 
This makes sense since, although we have computed A>z-yj at NLO, we should not expect 
improvement in the errors in the 0-jet bin without first calculating the NNLO corrections. 

In Figs. |] and || we present the dependence of the exclusive cross sections on the jet 
veto, for the 0- and 1-jet bins at NLO. We vary the scales using two different techniques. 
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Firstly we use the more traditional approach of choosing a central scale (in this case Mz) 
and varying it by a factor of two in each direction. Secondly we vary the scales using the 
approach discussed above (hereafter referred to as the ST method), in which the 0- and 
1-jet errors are defined as follows, 

(Af 7 LO ) 2 =(A^ 7 °) 2 + (A|« 7 ,) 2 , (5.4) 

(A^°) 2 =(A™,) 2 +(A|i^) 2 - (5-5) 

We calculate the uncertainties thus obtained for a range of veto scales. Note that, for the 
1-jet bin, we always require at least one jet with pt > 30 GeV and then vary the veto 
parameter for the second jet. It is clear from the figures that the ST method provides a 
more realistic measure of theoretical uncertainty than simply using the usual method of 
scale variation. The ST method has the pleasing feature of reproducing the inclusive error 
in the large-veto limit and, in addition, never results in a value of the jet veto for which 
the uncertainty vanishes. For example, this is the case for the Z7+jet calculation with 
pjp > 15 GeV (Fig. ||, right) when using a jet veto of around 30 GeV and the usual method 
of scale variation. 

Both figures also illustrate that the Z7+jet predictions have a larger ST scale variation 
than Zj. This can be explained by considering the differences between the Born production 
mechanisms in both. For the 0-jet bin the Born production process is purely electroweak 
and as a result the scale dependence is minimal, resulting from the factorization scale used 
in the PDFs. Clearly the dependence of the 1-jet bin on as occurs naturally at LO, thus 
yielding a much stronger dependence on (Xs(fJ-R) than the 0-jet bin. 

5.3 Z7 and Z7+ jet : Cross section summary 

To conclude this section we present the NLO predictions for the 0- and 1- jet exclusive 
and inclusive cross sections. We use the ST method to estimate scale uncertainty and 
also include uncertainties due to the PDFs and fragmentation contributions. The PDF 
uncertainties are obtained by using the 68% confidence level sets of CT10 [33]. In order to 
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Inclusive NLO (e+e~7 + X) [pb] 


pi > 15 GeV 
pi > 60 GeV 


1.21+°;^(scale)±3.5%(PDF)±0.4%(frag) 
0.0545+^J(scale)±3.7%(PDF)±1.2%(frag) 




Exclusive NLO (e + e~7 + no jets) [pb] 


pi > 15 GeV 
pi > 60 GeV 


1.02+}-^(scale)±3.8%(PDF)±0.5%(frag) 
0.0357+^ 4 ^ o (scale)±4.4%(PDF)±1.8%(frag) 



Table 2: NLO (inclusive and exclusive) predictions for the Z7 cross sections at LHC with ^fs = 7 
TeV. The N LO exclusive cross section, containing no identified jets, is defined using our usual jet 
cuts in Eq. (5.1). 



estimate the uncertainty arising from the fragmentation contributions we re-calculate the 
cross section using the fragmentation functions of Ref. [30] and compare the results with 
our default fragmentation set. This provides a crude estimate of the uncertainty arising 
from the modelling of the non-perturbative pieces of the fragmentation functions. These 
results are collected in Tables [2] and [3|. Note that imposing a jet-lepton separation (c.f. 
Eq. (|5.1|)) means that the inclusive Z7 cross sections presented in Table || depend on the 
jet definition. 

From these tables we can read off our theoretical predictions for the e + e~7 cross 
sections in the 0- and 1-jet bins, using the fiducial cuts employed in the ATLAS study [15]. 
Converting from relative to absolute uncertainties and adding them linearly we have: 

- jet : a NLO (pl > 15 GeV) = 1.02 ± 0.06 pb , a NLO (pl > 60 GeV) = 35.7±|;| fb 
1-jet : a NLO {pl > 15 GeV) = 187±|? fb , a NLO { P 7 r > 60 GeV) = 15.3±|;| fb 



In the 0-jet bin we see, by comparing with the predictions quoted in Eq. (52), that our 
revised results are similar to those presented in Ref. [15], with some of the difference at- 
tributable to the different choice of PDF set. However, the uncertainties are 50% larger due 
to the different treatment of the scale uncertainty. In the 1-jet bin, the combined theoreti- 
cal uncertainty is very large, indicating that a comparison with the theoretical prediction is 
of questionable value. In contrast the uncertainty on the inclusive 1-jet prediction (upper 
rows of Table [|) is much smaller, at the level of 10%, so that a much more meaningful 
comparison could be made. In passing, we note that our predictions for the 1-jet inclusive 
cross section are compatible with the difference between the inclusive and exclusive 0-jet 
bin Z7 results presented in Ref. [15]. 

A summary of our parton- level predictions, together with the ATLAS results from 
Ref. [15], is shown in Fig. ^. Note that all predictions are accurate to NLO, except for the 
2-jet results that are purely LO and therefore identical in the exclusive and inclusive cases. 
The uncertainty on the LO results corresponds only to scale variation, which is already 
considerable. 
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Inclusive 1-jet NLO (e+e"7 + jet + X) [fb] 


pi > 15 GeV 
pi > 60 GeV 


252+^g(scale)±2.0%(PDF)±0.6%(frag) 
24.6+JJ^(scale)±2.5%(PDF)±1.9%(£rag) 




Exclusive 1-jet NLO (e + e"7 + jet) [fb] 


pi > 15 GeV 
pi > 60 GeV 


188+^;°5(scale)±2.4%(PDF)±0.8%(frag) 
15.3+^;g(scale)±2.9%(PDF)±3.1%(frag) 



Table 3: NLO (inclusive and exclusive) predictions for the Z7 + jet cross sections at LHC with 
y/s = 7 TeV. Jets are defined using our usual jet cuts in Eq. (5.1). 
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p}>i5GeV pp^e + e"y 
p^>60GeV LHC 7 TeV 

i ATLAS 1 .02 fb" 1 
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p^> 15 GeV pp -> e e y 
p^>60GeV LHC 7 TeV 



1 2 
Number of jets 



ATLAS 1 .02 fb 
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Number of jets 



Figure 6: A summary of e + e~j cross sections obtained using MCFM. The left panel shows the 
exclusive predictions, i.e. for a specific number of jets and the right panel depicts the inclusive 
predictions, i.e. for at least the given number of jets. The ATLAS data, for the 0-jet bin in each 
case, are taken from Ref. fl5l. 



5.4 Photon pt spectrum in Zj and Z7+jet 

Next we consider the px spectrum of the photon produced in events containing ££7 and 
either zero or one jet. This is an important kinematic distribution since any deviation from 
the expected SM prediction may indicate the presence of anomalous couplings between the 
gauge bosons. In addition, searches in the missing £r+photon channel use the photon as 
a probe in order to search for the production of dark matter at colliders [35]. In both cases 
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Figure 7: The upper panel presents the photon pt spectrum using the low pt photon selection 
requirements at 7 TeV. Shown are the LO predictions (black dot-dashed), NLO exclusive prediction 
(blue-dashed) and NLO inclusive predictions (red-solid) . The left hand plot shows the photon pt 
for Z7, whilst the right hand plot is the spectrum from Z^-\-]et. The lower panel illustrates the 
inclusive (red) and exclusive (blue) if-factor for both predictions. 



accurate modelling of the background is essential in order to constrain, or observe, the new 
physics. 

Therefore in Fig. ^ we present the photon px spectrum for our low photon px selection 
requirement. We observe that in the region 15-200 GeV the inclusive A-factor is relatively 
stable in both the 0- and 1-jet bins. In the 0-jet case the K factor increases gently as the 
PT grows. However, for higher px values it approaches 2. The increasing A'-factor is hardly 
surprising since, as discussed in [5], at NLO one has large corrections to the inclusive rate 
from diagrams that include a gluon in the initial state, a high-pr jet and a relatively soft Z. 
Once this new kinematic regime is stifled by the application of a veto, the K factor reduces 
and becomes flatter, since the allowed kinematic region is more similar to the LO one. The 
inclusive A-factor for the 1-jet bin is smaller than its corresponding 0-jet counterpart. This 
is primarily because the presence of a jet in the Born topology allows more of the phase 
space for the Z and the photon to be explored at LO. In addition the presence of a gluon 
in the initial state at LO results in a more modest increase than the 0-jet bin since there is 
no large PDF enhancement of the real corrections. For the 1-jet bin the exclusive A-factor 
drops below 1 and falls significantly in the higher px region. We note that, for transverse 
momenta beyond the ranges that are presented here, our central scale choice of Mz may 
no longer be appropriate and an event-by-event scale such as <J M| + {p T ) 2 may be more 
reliable. 

It is interesting to consider the ratio of photon px spectra that could be constructed 
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Figure 8: The ratio of the photon transverse momentum spectra in charged lepton and neutrino 
final states for inclusive and Zj+jet production (left) and exclusive Z7 + 0-jet and Zj + 1-jet 
production (right). 



from a measurement of final states containing £ + £~j and l/P'y. This ratio may be useful for 
several reasons. Firstly the basic interpretation in the SM is that this ratio is sensitive to the 
contributions in which a photon is radiated from the final state leptons. Secondly, the ratio 
should suffer from fewer experimental ambiguities since one expects some cancellation of 
systematic errors. Thirdly, the ratio could be sensitive to models that modify the spectrum 
in only one of the channels. For instance, in dark matter scenarios only the photon spectrum 
in the missing Et + photon channel is modified. In contrast, in the case of anomalous 
couplings one would expect both px spectra to be altered in the same way so that the ratio 
is the same as in the SM. 

In Fig. H we present NLO predictions for such ratios in the presence of either or 1 
jets. In constructing these ratios we have considered a single lepton flavor but all three 
species of neutrinos. We observe that the 0— jet ratio has a strong peak in the low-p^ 
region associated with the radiation of softer photons from final state leptons. This is 
somewhat diminished in the 1— jet case, where there is an additional source of soft photons 
from events in which the Z boson almost balances with a hard jet. The tails in both cases 
tend to a constant ratio. We note that this constant is lower than might be expected from 
the relative branching ratios (~ 1/6), due to the difference in selection criteria for each 
process. For the invisible Z decay only the total neutrino momentum is subject to a pt cut 
and the rapidity is unconstrained, cuts that are less restrictive than those in the electron 
channel. 

5.5 Z77 phenomenology 

In this section we present some phenomenological studies for the final states t + £~^ and 
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1/^77 a t the LHC. Although the production cross section for Z7 is large the cost of radi- 
ating a further electroweak boson is severe, resulting in small rates at hadron colliders. As 
a result the observation of the Z77 cross section remains an experimental challenge. The 
process Z77 has received theoretical attention in the past, with NLO corrections calculated 
in Ref. [17]. However our study is the first to include the effects of photon fragmentation, 
allowing for the isolation of photons in a manner that is similar to that performed in 
experiments. We have checked our results against those presented in [17], using the Frix- 
ione isolation procedure described by Eq. ( [4,2|) , and find agreement within Monte Carlo 
uncertainties. 

Our aim in this section is to describe the phenomenology of the ^77 process, focussing 
primarily on the LHC operating at 8 TeV. Since the 8 TeV data set from 2012 alone could 
be around 20 fb _1 per experiment and the cross section for the e + e~77 process is about 
2 fb, analyses using this data set should have the best chance of observing this small SM 
process. The successful observation of this process at a hadron collider will at the very 
least instill further confidence in the ability of the LHC to identify rare processes with 
such small cross sections. In addition these rare SM processes may yield new insights into 
physics beyond the Standard Model. Observing a significantly different total rate than 
that which is predicted by the SM could be a sign of new physics. For this reason it is 
crucial to have predictions for the total rate accurate to NLO. 

We present cross sections for the production of l + £~~jj and z/1/77 using the following 
set of cuts, 

pi > 20 GeV , |7? 7 | < 2.5 , # 77 > 0.4 , 
p^ > 25 GeV , p e rf > 15 GeV , |%| < 2.5 , 
m u > 12 GeV , R h > 0.7 , Ef ss > 25 GeV , (5.6) 

where p^ ,s are the transverse momenta of the harder (h) and softer (s) leptons of the 
pair. These cuts are typical of those used in selection criteria at the LHC. For simplic- 
ity we maintain the same photon isolation requirements as in the previous section. For 
the common renormalization, factorization and fragmentation scale we choose [i = m^ 77 . 
Cross sections as a function of y/s are shown in Fig. [9|, together with the usual scale vari- 
ation. The total uncertainty on the cross section, including variations of both PDF and 
fragmentation sets, is rather small. For instance, the cross section for l + £~^ at 8 TeV is, 

a NLO (£ + £-jj) = 2.29 fb^5 (scale) ±4.1% (PDF) ± 0.5% (frag) . 

The K- factor across the range of operating energies considered is approximately 1.5 for the 
^ + ^~77 final state and 1.65 for i/z/77. The large iT-factors for these processes are reminis- 
cent of that obtained for Z7 [5] . In both cases the underlying Born process is qq- initiated 
and therefore the real corrections are significant, due to the large PDF enhancement from 
gluons in the initial state. 

We note that the cut we have used on m« is quite low, i.e we allow the lepton pair to 
be a long way from the Z pole. This benefits the analysis since there are a large number of 
£ + i?~77 events that contain at least one photon radiated from the leptons, thereby reducing 
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Figure 9: Dependence of the LO (blue) and NLO (red) inclusive cross section of pp — > £ + £~jj 
(left) and pp — > 3(vD)jj (right) processes on the renormalization/factorization scale, for a range of 
LHC operating energies. 



rri££. In Fig. 10 we present the invariant mass distribution of the two charged leptons using 
two values of the lepton-photon separation cut. We compare our usual choice presented 
in Eq. ( [5.6D , i.e. > 0.7, with the looser requirement > 0.4. From this figure it is 
clear that the invariant mass of the leptons is often away from the Z window. As expected, 
relaxing the isolation requirement dramatically enhances this region of phase space since 
the radiation peaks in the collinear region. These results suggest that the lepton-photon 
isolation should be as loose as is experimentally feasible in order to enhance the rate. 



6. Conclusions 

In this paper we have calculated NLO corrections to the production of Zj-\-jet and Z77 
at hadron colliders. We have included the full decays of the Z boson to leptons including, 
where appropriate, the radiation of photons from the Z decay products. We include frag- 
mentation contributions in order that photons can be isolated in a manner that is analogous 
to the current experimental procedure. Our results have been included in the Monte Carlo 
program MCFM that is publicly available. 

The results presented in this paper were obtained using analytic expressions for helic- 
ity amplitudes. In order to build these amplitudes the relevant Feynman diagrams were 
separated into gauge invariant subsets based on the identity of the fermion radiating the 
photon. This allowed an efficient recycling of earlier results for Zj and Zjj amplitudes 
into the necessary ingredients for the calculations at hand. The amplitudes in which a 
photon is radiated from a quark line are related to colour suppressed pieces of the qqZgg 
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Figure 10: NLO predictions for the invariant mass of the electron-positron pair in e + e~77 events at 
8 TeV. The dashed (blue) histogram shows the prediction using our usual electron-photon separation 
cut, i?£ 7 > 0.7 while the solid (red) prediction is for the looser cut, > 0.4. 

amplitudes originally calculated in Refs. [22, 36]. The remaining amplitudes, correspond- 
ing to photon radiation from the leptons, are presented here for the first time. They were 
obtained by extracting suitable QCD currents from Ref. [22] and then contracting them 
with the current for a Z boson decaying to leptons and one or two photons. 

We have studied the phenomenology associated with the production of a Z-boson in 
association with a photon and zero or one jets at NLO. The ATLAS collaboration has 
recently presented results for ££j cross sections [15], separating their results into bins 
classified by the number of jets present. Theoretical predictions for cross sections with a 
specified number of jets are subject to large uncertainties that can be underestimated by 
traditional estimates of the scale uncertainty [34]. We have shown that using the method 
of Ref. [34] indeed provides a more reasonable estimate of the theoretical uncertainty. We 
also consider both PDF uncertainties, by using the 68% confidence limit of CT10 [33], 
and fragmentation uncertainties, by comparing two independent fragmentation function 
calculations. We are thus able to provide theoretical predictions for the binned cross 
sections with uncertainties estimated using the best available information. We also studied 
the photon px spectrum in the various bins and presented ratios of pr spectra associated 
with charged and neutral leptonic decays of the Z-boson. 

Finally, we considered the much rarer Z77 process. Since the rate for this process is 
quite small it has so far not been observed at a hadron collider. In anticipation of a future 
measurement after the conclusion of the 8 TeV LHC data-taking, we have provided NLO 
cross sections and distributions for this process at that operating energy. 
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A. Amplitudes for Z7+jet production 

In this appendix we collect expressions for all the amplitudes required in the calculation of 
the Z7+jet process. 

A.l Tree amplitudes 

The decomposition of the tree level amplitude in terms of q- and l-type sub-amplitudes is 
given in Eq. (2.1). We begin by presenting the g-type tree-level amplitudes. As described 
in Sec. the expressions can be obtained from Ref. [22] by symmetrizing over the two 
colour orderings of the gluons in the amplitudes for — > qqggii. Explicitly we find, 

A q Hl<j~> 2 o » ^fl~J ^7 ' ' ^l) ~ 



(13) (14) (23) (24) (56) 



4°)(l+,2-,3+,4-,5r,6+) = -i (' ( 25 ) t 12 ] t 13 ] <5|(2 + 4)|3] + (12) (24) [16] (4|(1 + 3)|6] 



[14] [24] (56) t lu t 2U (13) (23) [56] t 134 t 2 34 

(4|(1 + 3)|6](5|(2 + 4)|3] (24) (25) [13] [16] 



[24] (13) ti 34 t234 [14] (23) ii 34 t234 



• (A.l) 



The tree-level ^-type amplitudes, obtained using the procedure described in Sec. are 
given by, 

a(°)(1+ o- 3+ 4+ 5r 6 + ) - ■ ^ 25 ^ 

\ L q i ^q i °g > ^7 j °£ > u £ 1 



(13) (23) (45) (46) ' 

-M-l- r i , . ■■ <2|(1 + 3)|6] 2 

t 456 (13) (23) [45] [46] 



4 0) (1^2r,3+ 4- 5J, 6 +) = i - ■ (A- 2 ) 



A. 2 Virtual amplitudes 

We again begin with the decomposition of the — > qqgjll one-loop amplitude into sub- 
amplitudes, 

A«(l 9 ,2 9 -,3 9 ,4 7 ,5,-,6,) = 2v / 2e 3 g^crT^ 

Q q (~Q q + ^Ir'Pz(s^)) A$0. q , 2 q , 3 g , 4 7 , 5,-, 6,) 

+Q q £ ( ~Q< + 2 V ln( v l + v R )Vz{s^) ) A%(l q , 2 f , 3 9 , 4 7 , 5,-, 6,) 

i=l ^ ' 

v e 

+Q q 7T^ V zMA q 1 lO- q , 2 q , 3 g , 4 7 , 5j, 6*) 
+Qe (-Q q + v q L R vl R V z (t^)) 4?i ) (l< ? , 2 5,3 9 ,4 7 ,5 £ -,6,) 

+Q£^ ± ^^z(t45 6 )^ 1 2 ) (l g , 2 ? , 3 S , 4 7 , 5,-, 6,)) , (A.3) 
2siyejy ' J 
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where the usual one-loop prefactor cr is defined as, 

1 r(i + e )r 2 (i-e) 



cr 



(4vr) 2 -^ r(l - 2e) 



(A.4) 



The g-type contributions are, 



A^l(l q ,2 q ,3 g ,4^,5 b 6i) = -iV c A 6 (l ? ,3,2 ? -,4) 



A 6 (l ? ,2 g -,3,4) +A 6 (l 9 ,2 9 -,4,3) 



"^g;2(^<?» ^<?» ^gi ^£1 ' 



AI 4 (l 9 ,2 9 -,3,4)+^. 4 (l ? ,2 ? -,4,3) 



2 q , 3 g , 4 7 , 5,-, 6,) = ^ 4 (1„ 2,-, 3, 4) + A&{l q , 2 q , 4, 3) 



(A.5) 



where below we list references for the primitive amplitudes from the Ref [22], 

A 6 (l g ,2,3 5 ,4) 
A 6 (l g ,2 ? -,3,4) 

^■6;4(^9' 4) 

^6*4 (Igi ^g, 3, 4) 
The ^-type contributions in Eq. 



Sec. 9, Eqs. (9.2)-(9.14), 
Sec. 10, Eqs. (10.1)-(10.21), 
Eq. (2.13), 
Eq. (2.13). 



are, 



A { 3{l q , 2g, 3 g , 4 7 , 5,-, 6,) = iV c ^ c (l g , 2g, 3 S , 4 7 , 5,-, 6,) + -^AfO-q, 2 q , 3 g , 4 7 , 5 



iV, 



-^£ ; 2(l(j) 2 g , 3 3 , 4 7 , 5^, 6^) — ^ (l g , 2 g , 3 5 , 4 7 , 5^ 



(A.6) 



where the sub-amplitudes have been divided into three contributions: leading color (lc), 
subleading color (si) and fermion-loop (fl). We now present analytic formulae for the 
three contributions appearing in Eqs. ( |A.6| ). We again follow the notation of Ref. [22] and 
decompose our one-loop amplitude into divergent (V) and finite (F) pieces, 



At = V e A? } + %F t 



(A.7) 



We will present results for the two independent helicity combinations, (1+, 27, 3+, 4+, 57, 6+) 
and (1+, 27, 3+, 4~ , 57, 6+). The other two gluon and photon helicity combinations, for the 
same quark and lepton helicities, are obtained from the following relations, 



.4^1+ 27,37,4+ 57,6^ 
At(lq,2q- ,37,47/ ,5j,6j 



A e (2± 17,3+47,67,5^ 



A e (2+ 17, 3+4+ 67, 5j 



For the case where the helicity of the lepton pair is flipped, the amplitudes are obtained by 
performing a 5 o 6 exchange with an extra minus sign due to the Z/'y* —> l + £^j current, 



5^6 



Ae(l q ,2q, 3 g ,4 7 ,5^ ,6 £ ) = -Af(l q ,2q, 3 g ,4 7 ,5^ ,6 £ ) 
Similarly, the amplitudes for the remaining helicity combinations can be obtained via, 
AeO-q ,2+,3 9 ,4 7 ,5^-,6£) = -Ae(lq, 27, 3 S , 4 7 , 5|, 6^)| lo2 . 
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All helicities share a common divergent factor, 
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For the helicity configuration (1+, 2g , 3+, 4+, 5^ , 6^) the finite remainder is, 



Fl c = - 



(25)' 



(13) (23) (45) (46) 



Ls_ 



-«13 -S23 



-^456 —^456 



(12) (25) (5|(4 + 6)|1] ^\=u£ 



(12) 2 (5|(2 + 3)|1](5|(4 + 6)|1 



L 



i - 



<456 



2 (13) (23) (45) (46) 



"-456 



(13) (23) (45) (46) 



^456 
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and for (1+ 2 q , 3+ 4 5- , 6 



(2|(1 + 3)|6](2|(4 + 5)|6] 
i 456 (13) (23) [45] [46] 



Ls 



si3 S23\ (2|1|6](2|(4 + 5)|6]) 



"^456 —^456 



-^23 



-<456 / 



(13) (23) [45] [46] i 456 



[12] 2 [16] z 



L-, 



-323 



"*456 / 



2 (13) (23) [45] [46] t\ m 
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Explicit formulae for the basis integrals appearing in this formulae can be found in Ap- 
pendix II of Ref. [22]. 

Next we consider the sub-leading colour contribution to Ai{V^ , 2^, 3+, 4+, 5J, 6^). Us- 
ing the same notation as before the divergent pieces are given by, 



Vf 



+ 



2e \— £456 



+ 



(A.11) 



The finite part for the (1+, 2^ , 3+, 4+, 5^ , 6^~) configuration is. 



TlSl 



(25)' 



(13) (23) (45) (46)' 



-Ls_ 



S\1 -«13 



-t456 — *456 



+ 



(12) 2 (35) s 



(13) 6 (23) (45) (46) 



-Ls_ 



-syi -S23 



-t456 — *456 



2ai 3 <15) (25) - (2|3|1](15) 



2L °(^f) 



[13] 2 (23)(15) 



2 Li 



(13) 2 (45) (46) 



S23 



2 (13) (45) (46) 



'23 



(2|1|3](15) (35) 



+ 



(2|1|3](53) <5|(4 + 6)|3] 



(13) z (45)(46) si2 (13) (45) (46) 

(5|(1 + 2)|3] ([13] (5|(4 + 6)|2] + [23] (5|(4 + 6)|1]) 
2i 456 [12] [23] (13) (45) (46) 



-^456 



'12 
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while the corresponding contribution for (1+ 2= , 3+ 4 , 5^ , 6^") is, 



, (2|(1 + 3)|6](2|(4 + 5)|6] 
1 t456 (13) (23) [45] [46] 1 

(12) 2 (3|(1 + 2)|6](3|(4 + 5)|6] 
t 456 (13) 3 ( 23) [45] [46] 



q ' q ' "g ' 7 ' i 

Sl2 Sl3 



Ls_i 



£456 — *456 

-Sl2 -S23 



"^456 — *456 



(1|(2 + 3)|6](2 Sl3 (2|(4 + 5)|6] - (2|3|1](1|(4 + 5)|6]) 



t 456 (13)" [45] [46] 

[13] 2 (23)(1|(2 + 3)|6](1|(4 + 5)|6] L i (^ff) 
2t 456 (13) [45] [46] s| 3 

(2|1|3](3|(1 + 2)|6](1|(4 + 5)|6] L o (^) 



S23 



(A.13) 



t 456 (13) 2 [45] [46] 
[36] ([13] [62] + [23] [61]) 



«12 



+ 



(2|1|3](63) (3|(4 + 5)|6] 
(13) [45] [46] 



( ~*456 ) 



'12 



2 [12] [23] (13) [45] [46] ' 

Lastly we consider the fermion loop contribution, which does not contain a divergent piece 
(V/ = 0). The finite part for the (1+, 2^, 3+, 4+, 5J , 6+) configuration is 



fl = [13] (25) (5|(4 + 6)|3] 
£ " (45) (46) 



Li 



-Sl2 
-<456 



1 



t 2 



12t 456 77i| 
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and the result for the the other helicity configuration (1+, 2^ , 3+, 4 7 , 5^ , 6^) is, 



rfl = [13] [36] (2|(4 + 5)|6] 
1 " [45] [46] 



Li 



-S12 \ 
"<456 / 



1 



12t 456 m/ 
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Here the top quark loop contribution is expanded in powers of 1/ m 2 and the terms of order 
1/mf are dropped. 

A. 3 Real emission amplitudes 

Next we consider the real corrections to the Z7+jet process. We use the same notation as 
in the previous section to designate q- and i- type diagrams. 

A. 3.1 Amplitudes for — > qqggjU 

We begin by considering the amplitudes containing only one quark line, i.e. the process 
— > qqgg^tl. The decomposition of this amplitude is given by, 



A(°>(1„2 9 -,3 9 ,4 S ,5 7 ,6*7*) = 2V2e 3 g 2 s £ (T a *T a % 

{3,4} 



■l 2 



Q, (~Q q + <,rA,rVz{s q7 )) 4°) (1,, 2g, 3g, 4g, 5 7 , 6,-, 7() 
+Qe (~Q q + v% R vi iR Vz(t 5(i 7)) Af ] (1, , 2 5 , 3 g , 4 g , 5 7 , 6,-, 7,) j . (A.16) 
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All that has changed with respect to the tree-level decomposition is the presence of the 
additional gluon, which manifests itself in the form of the amplitudes and the colour struc- 
ture. 

The helicity amplitudes for diagrams in which the photon is emitted from the quark 
line are obtained from the amplitudes for the e + e~ — > qqggg process presented in Ref. [36], 
where one gluon is replaced by a photon. 



4°)(l 3 ,2,-,3 g ,4 g ,5 T ,6^7^) = — 

S67 



^-(1<?> Ag, 5g, 2 q ) + A(l q , 3g, 5g,4g, 2 q ) 



+ A(l q , 5g, 3g, Ag, 2 q ) 



(A.17) 



where on the right-hand side, the 6, 7 labels for lepton pair have been suppressed. The 
analytic expressions for A{l q , 3 g , 4 g , 5 S , 2 q ) are presented in Appendix A, Eqs. (A43)-(A49) 
of Ref. [36]. The helicity amplitudes for diagrams where the photon is emitted from the 
lepton line are given by 

4 0) (l+,2r,3+,4+,5+,6r,7+) = i ^ J^J^ (5g) (57) , (A.18) 
.(0) a+ + 4+ 7+) _ (2|(5 + 6)|7] 2 

.(0) f + 4 _ + 7+) _ (6|(5 + 7)|1] 2 

JK t {i q , Z q , 6 g , 4 g , t> 7 , b ? , ( e ) - [i3] [34] [42] (5g) (57) , [A. ZD) 

(o) [i 1 ^]^ 

V (1+ 2r,3J,4J, 57,67, 7^) = -i [l3] [M] m [5g] [5?] , (A.21) 

(0) f + + + + , _ i f (4|l|3](26) (4|(l + 3)(5 + 7)|6) 

A e [l q , 2, , 6 g , 4 g , 5 7 , b, , 7, ) - s34t5g7 (56) (57) I (13) 

(4|2|3](6|(4 + 2)|3](6|(5 + 7)|1] (4[(1 + 3)(5 + 7)|6)(6|(4 + 2)|3] ] 

[42}t 2U (13) [42] / ' { ■ ' 

(o )fl+ 2 _ 3+ fi _ ?+ , _ i f (4|l|3](4|(l+3)|7](2|(5 + 6)|7] 

1 { q ' q SMt567 [56] [57] { (is) tl34 

(4|2|3][3|(4 + 2)(5 + 6)|7] [17] [3|(4 + 2)(5 + 6)|7](4|(1 + 3)|7] ) 

[42]t 234 (13) [42] /' ^ 

A (0) (1+ 2 - 3- 4+ 5+ 6 - 7+1 _ i f [14] 2 (26)(3|(l + 4)(5 + 7)|6) 

^ ^ , Z- q , 3 ff , 4 g , 5 7 , b, , 7 e ) - ss4t567 (56) (57) I [ig] ^ 

(32) 2 (6|(3 + 2)|4](6|(5 + 7)|1] [14] (32) (26) (6|(5 + 7)|1] | 
+ (42) t 234 + [13] (42) J ' (A - M) 

(0 ) ri+ 2 _ 3- 4+ 6 _ 7+ , _ t \ [14] 2 (3|(1 + 4)|7](2|(5 + 6)|7] 

(32) 2 [17][4|(3 + 2)(5 + 6)|7] [14] (32) [71] (2\(5 + 6)\7} \ 
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A.3.2 Amplitudes for -)• qqQQ-iU 

Next we consider the processes that contain two quark lines. The decomposition of the 
amplitude is, 



A^(l g , 2g, 3q, 4q, 5 7 , 6,-, 7,) = 2V2e 3 g s T^T^S^d QQ 
{-Q q + v% R vt iR Vz(s67)) Af\l q , 2g, 3 Q , 
+Qe {~Q q + v q L)R vi A V z {h G7 )) Af\l q , 2 q , 3 Q , 4q, 5 7 , 6 f , 7^ 
+ ^2,^1,^(^67)) 4 0) (3q, 4q, 1„ 2 9 -, 5 7 , 6,-, 7,) 



(A.26) 



+Q< ( "Qq + v^ R v% R r z (t 567 ) ) ^ 0) (3 Q , 4 , 1,, 2,, 5 7 , 6,-, 7,) 



2? 4 <; 



The four quark amplitudes are obtained from Ref. [36] in a similar fashion as described 
above. The helicity amplitudes for diagrams where the photon is emitted from the quark 
line are given by, 



4°)(l 9 ,2 9 -,3 Q ,4 ,5 7 ,6,-,7,) = — 



$67 



Qq [Al(l q ,2q,3 Q ,4:Q,5 g ) + A 2 (l q ,2q,3 Q ,4Q,5g)) 



+QQ [A 3 (l q , 2 q , 3 Q , 4q, 5 g ) + A 4 (l q , 2 q , 3 Q ,A Q , 5 g )) 



(A.27) 



The analytic expressions for Ai(l q , 2 q , 3q,4q, 5 g ) are presented in Appendix A, Eqs. (A50)- 
(A59) of Ref. [36]. For diagrams where the photon is emitted from the lepton line the 
helicity amplitudes are, 



a(°)(-\+ o- 2+ 4~ 5+ 6- 7 



— I 



9 ' S34 t 567 (56) (57) 

[13](26)(4|(l + 3)(5 + 7)|6) , (42) (6|(2 + 4)|3](6|(5 + 7)|1] 



1 *134 

a(P)(1+ 2- 3 + 4~ 5~ 6~ 7 



+ 



*234 
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•534*567 [56] [57] 



[13](4|(1 + 3)|7](2|(5 + 6)|7] , (42) [71] [3| (2 + 4)(5 + 6)|7] 



+ 



H34 



*234 
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B. Amplitudes for Z77 production 



The decomposition of the tree level — > qqjj££ amplitudes in terms of qq-,q£- and ii- 
type sub-amplitudes is given in Eq. fl3.1| ). These sub-amplitudes can be obtained from the 
results of the previous section using the following relations, 



A$(l q , 2 q , 3 7 , 4 7 , 5 b 6 e ) = A g °\l q , 2 q , 3 g , 4 7 ,5 £ -, 6 e ) 
A q e 0-q, 2 q , 3 7 , 4 7 , 5^, 6e) = A\ 2 q , 3 g , 4 7 , 5^, 6^) 
A^(^-q, 2 q , 3 7 , 4 7 , 5^, 6^) = .4.^(6^, 5^, 3 7 , 4 7 , 2 q , l q ), 



(B.l) 
(B.2) 
(B.3) 



- 26 - 



The amplitudes A q °^ and Af are denned in Eqs.(A.l) and ( |A.2| ) respectively. Note that 
Eq. flB.2j ) represents the tree-level amplitude for the emission of photon 4 from the lepton 
line. 



Due to the simple colour structure, the extension of Eq. ( p.l| ) to the one-loop case is 
simple, 

AW (i g , 2g -, 3 7 , 4 7 , b b 6,) = c r g 2 s (n c - ±?j A^(l q , 2 q , 3 7 , 4 7 , 5* 6,) {^y -> 4y } ■ 

(B.4) 

Analogous to the Z7+jet process, the helicity amplitudes for the gg-type diagrams can be 
obtained from the amplitudes for the e + e~ —> qqgg process presented in Ref. [22]. In par- 
ticular we only need the subleading color primitive amplitude Ae(l q , 2 q , 3,4) symmetrized 
over the two gluons. The q£-type amplitudes can be obtained in a similar fashion from the 
subleading color contribution to At(l^, 27, 3+, 4+, 57, 6^T) amplitude, which we presented 
in Appendix |A]. Specifically in terms of the amplitudes presented in [22] and the previous 
section we have, 



*^qq 0-q> 2q, 3 7 , 4 7 , 5^, 1 



A 6 (l„2 ? ,3,4) + A 6 (l„2 ?) 4,3) 



^ ( 1 g> 2 g ) 3 7 ,4 7 ,5^,6^) = (1 9 , 2 f , 3 S , 4 7 , h h 6^) . 



(B.5) 
(B.6) 



The £^-type contribution to the — > qqjj££ one-loop amplitude simply consists of a 
vertex correction and is given by 



■AitO-Q) 3 7 , 4 7 , 5^, 6i) 



-S12 



-1-1-4 
e 2 2e 



4?(l g ,2 9 -,3 7 ,4 7 ,5,-,6,). (B.7) 



Finally we consider the real corrections in which an additional gluon is radiated relative 
to the Born process. In terms of our usual decomposition the —> qqg^ll amplitude is 
given by, 



A<®(l q , 2- q , 3 g , 4 7 , 5 7 , 6 E , 7 t ) = ie^T^ 
x Q 2 q (~Q q + v% R vi R V Z (s 67 )) A$(l q , 2 q , 3 g , 4 7 , 5 7 , 

%, 7i) 



+QeQ q [~Q q + vl R vl R V z (t467)) A^ (l q , 2 q , 3 g , 4 7 , 5 7 , 6* 7. 
+QiQ q [-Q q + v q L R vi iR Vz{tmi)) A$(l q , 2 q , 3 g , 4 7 , 5 7 , 6 b 7. 
+Qj \ -Qq + v\ R vl !R V z (h567)) A% ' (l q , 2 q , 3 9 , 4 7 , 5 7 , 6^, 7 e ) 



(B.8) 

The sub-amplitudes can be constructed as follows, 

(1„ 2 q , 3 g , 4 7 , 5 7 , 6* 7 t ) = A q 0) (l q , 2 q , 3 g , 4 S , 5 7 , 6* 7<) + A q 0) (1„ 2 5 , 3 S , 5 7 , 4 9 , 6* 7*) 



•Arf (1<2' 2?) 3pj 4 7 , 5 7 , 6^, 7^) — .A^ ; (lq, 2 q , 3 g , 4 ff , 5 7 , 6^, 7^) + .A^ ; (1 9 , 2 9 , 4 9 , 3 9 , 5 7 , 6^, 7^) 



,(o) 



(o) 



V 
,(o) 



•^ti 2g, 3 S , 4 7 , 5 7 , 6^, 7^) — Aq£ (7^, 6^, 5 S , 3 7 , 4 7 , 2 q , 1 



(B.9) 



where A^ (l q , 2 q , 3 g , 4 S , 5 7 , 6^, 7<?) and 2,, 3 S , 4 ff , 5 7 , 6^, 7^) are given in Eqs. ( |A.17| ) 

and ( |A.18|) - (|A.25| ) respectively. 
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